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We consider the phase diagram of QCD at very high baryon density and at zero temperature 
in the presence of a strong magnetic field. The state of matter at such high densities and low 
temperatures is believed to be a phase known as the color-flavor locked phase which breaks color 
and electromagnetic gauge invariance leaving a linear combination of them, denoted as U(l) e fn 
unbroken. Of the 9 quarks (three flavors and three colors), five are neutral under this unbroken 
generator and four are oppositely charged (two with a charge of +1 and two with —1). In the presence 
of a magnetic field corresponding to U(l) e fn however, the properties of the condensate changes and 
a new phase known as the magnetic color flavor locked (MCFL)phase is realized. This phase breaks 
an approximate SU(2)c x SU(2)l x SU(2)r x U(1)b x U(1)a symmetry of the Lagrangian to 
SU(2)c+l+r giving rise to 6 Goldstone modes 5 of which are pseudo Goldstone modes. These 
Goldstone modes are composed of excitations that correspond to both neutral quarks and charged 
quarks. Hence it is natural to expect that the propagators of these Goldstone modes get affected in 
the presence of a magnetic field and their speed becomes considerably anisotropic. Although this 
anisotropy is self-evident from symmetry arguments, it has not been quantified yet. We calculate 
this anisotropy in the speed of the Goldstone modes using an NJL model type of interaction between 
the quarks and comment on the impact of such anisotropic modes on transport properties of the 
MCFL phase. 


I. INTRODUCTION 

The low energy properties of any state of matter is primarily dictated by its lowest energy excitations. In case of a 
Bose condensate these excitations are usually the gapless Goldstone modes which arise due to the spontaneous break¬ 
ing of continuous symmetries by the condensate. This paper explores the effect of a magnetic field on the propagation 
of the Goldstone modes that arise due to the spontaneous breaking of color-flavor symmetries in a condensed phase of 
matter known as the color-flavor locked (CFL) phase in the presence of a magnetic field. Before going into the details 
of what color-flavor locking is we can expect to predict some of the consequences of having a magnetic field using 
symmetry arguments alone. The presence of a magnetic held breaks rotational symmetry explicitly and the most 
natural consequence of it on the propagator of the Goldstone modes can be a difference in the speed of the modes in 
a direction perpendicular to the Held compared to a direction parallel to the Held. Although the effective theory of 
the Goldstone modes of the CFL phase in the presence of a magnetic field has been discussed in some detail in [1], 
the effect of a magnetic field on the speed has been commented on at a qualitative level. We attempt to quantify the 
effect of the magnetic field on the neutral Goldstone mode propagators of color-flavor locked phase in this paper. 

To put things in context, the phase diagram of QCD at finite baryon densities garners a lot of attention due 
to its implications for the physics of compact astropliysical objects, terrestrial nuclear experiments and heavy ion 
collisions. Although considerable progress has been made in exploring the phase diagram, some parts of it still remain 
inaccessible to rigorous analytical or numerical treatment. The region of the phase diagram that corresponds to 
moderate to very high baryon densities at low temperature is particularly relevant for the physics of compact stars, 
the reason being that these stars are cold and composed mostly of baryons. However, the regime with moderate 
baryon density and low temperature is difficult to analyse due to the following reasons. The strength of the QCD 
coupling constant at low energies being large renders perturbative calculations from first principles futile necessitating 
the use of numerical calculations on lattice. But lattice simulations at finite baryon density are not useful either as 
the lattice algorithms involving important sampling break down due to the sign problem[2-6]. However there is a 
region of the phase diagram, at very high baryon densities and low temperatures, that can be dealt with analytically. 
This can be explained as follows. At high density and low temperature quarks form Fermi spheres. The physics, 
in this regime of the phase diagram, is mostly dictated by the quarks at the Fermi surface. As the fermions at the 
Fermi surface are highly energetic, they are weakly coupled making rigorous perturbative calculations justified. This 
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TABLE I: The t/(l) e m charge of the quarks 


regime also happens to be relevant for the physics of the core of neutron stars as the density of baryons at the core is 
expected to reach asymptotically high values. The high density regime has been explored in great detail in the past 
two decades. In this regime matter was predicted to exist in a colorsuperconducting phase known as the ‘color-flavor 
locked’ phase using both model calculations [7-12] and weak coupling calculations in QCD [13-17], [18-20] and 
[21-23]. The term ‘color superconductivity’ corresponds to the Bose condensation of color-charged Cooper pairs of 
two quarks at the Fermi surface. As mentioned before at high densities and low temperatures quarks form Fermi 
spheres with a radius that depends on the baryon chemical potential. In such a scenario, if there is an attractive 
color-flavor and angular momentum channel available at the Fermi surface giving rise to a BCS instability, two quarks 
at the Fermi surface can form a Cooper pair lowering their energies. The Cooper pairs being Bosonic excitations, 
form a Bose-Einstein condensate which is not color neutral. The gluons acquire a Meissner mass in such a condensate 
and get screened. Many such color superconducting phases with different quark pairings have been proposed for a 
range of baryon densities. At asymptotically high densities the energetically favoured color superconducting phase 
was found to be the color-flavor locked phase. The condensate is not color neutral as expected and breaks both color 
and U( l)em- However, there is a linear combination of the 8 th generator of the color SU C ( 3) and the U( l) em that 
remains unbroken by the condensate. This linear combination is referred to as the rotated electromagnetism and 
denoted by U(l) e fn- A magnetic field corresponding to this rotated electromagnetism (B) does not get screened by 
the condensate. The scenario at this point diverges markedly from the U( l) em superconductivity in terrestrial metals 
which shields magnetic fields completely or accomodates magnetic fields only inside quantized vortices. This opens 
up exciting new possibilities as far as the properties of CFL phase in the presence of a magnetic field is concerned. 
But large magnetic fields in CFL phase are interesting not for theoretical reasons alone. There is usually a very 
strong magnetic held associated with a compact star. For a neutron star the typical fields on the surface are such that 
B ~ 10 12 Gauss and for a magnetar the fields could be as large as 10 16 Gauss[24-26]. For gravitationally bound stars 
an upper bound on the strength of the magnetic held can be found by comparing the energy of the gravitaional held 
and that of the magnetic held of the star. This bound is around 10 18 Gauss. If we consider self bound quark stars 
however, this upper limit can be even higher [27]. Also, higher magnetic helds are expected to exist inside [/(l) e m 
charged gluonic vortices [28, 29] which arise due to chromomagnetic instability in dense matter[30-32] in the presence 
of a magnetic held. Typical magnetic held strength in the core of magnetars may well reach the order of the square 
of the CFL gap or higher and can be expected to alter the nature of the diquark pairing of the condensate giving rise 
to a different color superconducting phase known as the magnetic CFL or MCFL phase. 


The symmetry properties of the magnetic CFL phase are slightly different from that of the CFL phase. To un¬ 
derstand what these differences are let us hrst look at the symmetry breaking patterns of the CFL phase in the 
absence of a magnetic field. The standard model Lagrangian with a baryon chemical potential has approximate 
SU( 3) c x SU(3) L x SU(3) r x U(1) b x U(1)a symmetry well below the weak scale. The condensation of CFL phase 
breaks it to SU(3)c+l+r giving Meissner mass to seven of the eight gluons. A linear combination of the 8 th gluonic 
generator and the generator of U{ l) em also becomes massive where as the orthogonal linear combination of the two 
remains massless (U(l) e m)- The CFL phase has a total of 10 Goldstone modes including 7r ± , ir°, K°,K °, A' + , K~, rj, rj 
mesons in the massless quark limit. In addition there is a Goldstone mode due to the breaking of U(1) B which we 
denote by </>. In order to predict the observable consequences, if any, of the existence of a CFL phase in the core 
of a neutron star it is of utmost importance to analyse the effective theory of these Goldstone modes. The low 
energy effective theories of these Goldstone modes were discussed in great detail in [33-38] more than a decade 
ago. The effect of a XJ{l) e m magnetic field B on these Goldstone modes however remains to be explored. But 
in order to discuss the effects of B on these Goldstone modes adequately we need to also take into account how 
the pairing of the condensate and consequently the symmetry properties get affected by this magnetic held. In 
the presence of B, the aforementioned SU(3) C x SU(3)l x SU(3)r x U(1)b x U{\)a symmetry of the standard 
model Lagrangian gets explicitly broken down to SU(2) C x SU{2) L x SU(2) R x U(1)b x U(1 )a symmetry as four 
of the nine quarks are electrically charged under this U(l) e m- The charges of the quarks are shown in Table. I. 
When the magnetic held is small compared to the CFL gap the effect of the magnetic held on the structure of the 
gap is negligible and the form of the condensate is given by that of the CFL pairing. However, when the mag¬ 
netic held becomes comparable to the CFL gap, the form of the quark-pairing changes and the condensate breaks 
SU(2) C x SU(2 ) l x SU{2 ) r x U{l) B x U( 1) A down to SU(2) R+L+C [27, 39, 40]. The gap in the MCFL phase was 
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calculated for both magnetic fields that are extremely large (/z 2 = eB , ~ 10 19 Gauss for /z ~ 500MeV) [27, 39, 40] 
and for smaller magnetic fields (10 19 Gauss > eB > 10 18 Gauss) [27]. It was also shown that if the magnetic field 
is lowered further the gap structure of the magnetic CFL phase becomes identical to that of CFL. The lower limit 
Four of the Goldstone modes of the CFL phase are massive in the MCFL phase due to the explicit breaking of some 
of the color-flavor symmetry of the QCD Lagrangian. These are the modes that are charged under the U(\)em and 
are basically the charged pions and kaons. The MCFL phase retains the other 6 Goldstone modes of the CFL phase 
listed above. The effective theory of these Goldstone modes was discussed to some extent in [1]. These six modes are 
neutral under the rotated U(l) e fn and are made up of either two neutral quark excitations or two oppositely charged 
quark excitaions so as to keep the combined excitaion neutral. Small magnetic fields do not affect the Goldstone 
mode propagation appreciably as the modes are neutral under this U(l) e m- However a large magnetic field B should 
be able to resolve the internal structure of the modes and be able to change their properties. One such effect is 
encapsulated in the speeds of the Goldstone modes and we can quantify this effect by looking at the underlying 
theory from which the effective theory of the Goldstone modes is derived. In the presence of a magnetic field B 
the quarks are basically stuck in Landau levels which not only alters the gap structure of the MCFL phase from 
the CFL phase but also causes the speeds of the neutral Goldstone modes to be different in directions parallel to 
the magnetic field compared to that perpendicular to it. The six neutral Goldstone modes are of interest to us as 
being neutral they are massless in the presence of a magnetic field and easy to excite at temperatures much below 
the gap. This means that the low energy properties of the MCFL phase are dictated by these neutral Goldstone modes. 

The paper is organized as follows. In section Sec. II we analyse magnetic fields as high as eB ~ 2/z 2 and in 
section Sec. Ill we look at more realistic smaller magnetic fields followed by a concluding section. Although eB ~ 2/z 2 
may be too high a magnetic field to be present at the core of a compact star, it is not ruled out and may well be 
within the range of realistic magnetic fields in quark stars. Also, a mechanism was conjectured [28, 29]by which high 
magnetic fields can be found to exist inside gluonic vortices in CFL phase. The mechanism is as follows. Under 
the rotated U(l) e m, some of the gluons are charged. This causes the dispersion relation of one of the polarizations 
of these charged gluons to become tachyonic for high magnetic fields creating a chromomagnetic instability. As a 
consequence charged gluons condense in a vortex state with magnetic flux tubes inside, which correspond to stronger 
magnetic fields than the magnetic field present outside the vortices. This phenomenon is called antiscreening of the 
magnetic field and can give rise to magnetic fields in the core of magnetars that are stronger than expected otherwise. 
Despite the magnetic field being too large, it is worth analysing such fields as the computation of the propagator for 
the neutral Goldstone modes is simplified and enlightening for the purpose of demonstrating how a similar calculation 
will work for smaller magnetic fields. 


II. GOLDSTONE MODE PROPAGATOR AT VERY HIGH MAGNETIC FIELDS 


Before we compute the propagators for the Goldstone modes, the energy scales of interest need to be specified. 
Here, we are interested in the neutral gapless excitations about the MCFL condensate with momentum much smaller 
campared to the gap which is much smaller than the chemical potential in the presence of a magnetic field which is 
of the order of the square of the chemical potential. This means that the quarks occupy only the lowest landau level. 
We start with a Lagrangian with 3 colors and 3 flavors of quarks interacting via an NJL type four-Fermi interaction 
in the color-flavor locking channel in the presence of a U(l) e fn. gauge field with a baryon chemical potential /z 
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C = ip(i$ + eQA + /z7o)V> + ^ ^(ipP v ip c )(ipcP v ip) (1) 

ri—l 


where z/> ia is the quark field with a being the flavor index and i being the color index. Also the field ip c is the charge 
conjugate of the field ip given by ip c = Cip T where C = i'y 2 j°. A M is the U(l) e fn gauge field. For simplicity we restrict 
ourselves to the anti-symmetric color-flavor locked channel. Hence, (P v ) ia _j b = i^btij^abr)- Q is the rotated U(l) e fn 
charge of the quarks and is given by (+1)U + + (—l)f2_ + (0)U o where U + ,U_ and f2 0 are projectors of positively 
charged quarks, negatively charged quarks and neutral quarks. The operators 0±,Uo can be expressed as follows in 
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a notation that will be explained in the few lines that follow, 


/0 

0 

0 

0 

0 

0 

0 

0 

°\ 


/o 

0 

0 

0 

0 

0 

0 

0 

o\ 


(1 

0 

0 

0 

0 

0 

0 

0 

0\ 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

1 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

, n . = 

0 

0 

0 

0 

0 

0 

0 

0 

0 

, fi 0 = 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

1 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

Vo 

0 

0 

0 

0 

0 

0 

0 

0/ 


\0 

0 

0 

0 

0 

0 

0 

0 

0^ 


\0 

0 

0 

0 

0 

0 

0 

0 

1/ 


To clarify our notation, for i = ±,0, is written as a 9 x 9 matrix above, where, 3x3 combinations of 3 

flavors(a = 1,2,3) and 3 colors (i = 1,2,3) are denoted by the 9 rows of the matrix and 3x3 combinations of 3 
flavors(& = 1,2,3) and 3 colors (j = 1,2,3) are denoted by the 9 columns of the matrix. G is the coupling constant 
with dimension of length 2 . As we are interested in the bosonic excitations about the condensate, we need to introduce 
auxiliary fields in the action. Remembering that we are considering pairing only in the anti-symmetric channel in 
color and flavor, we introduce three auxiliary fields Ai, A 2 and A 3 and rewrite the action as 


The matrix ;S' quark is given by 


S = — 


0-1 

quark 


A 2 + A 2 


G 
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where A represents the auxiliary fields in our 9x9 matrix notation as 
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( 5 ) 


We introduce an external magnetic field at this point in the problem by assigning = (0, 0, Bx, 0). The idea is 
to first minimize the action of Eq. 3 with respect to the magnitude of the auxiliary fields which gives us three gap 
equations from which we can extract an estimate of the magnitudes of the three gaps. Then we expand the action of 
Eq. 3 about this minimum, to quadratic order in the fluctuations to obtain the leading order terms of the effective 
action for the Goldstone modes. For simplicity we choose to work in a regime where the three gaps are equal, i.e, 
Ai = A 2 = A 3 = A. But in order to perform a consistent calculation we need to figure out where in the QCD phase 
diagram (/i vs B) this regime is or at what value of the magnetic field the three gaps equal each other. But before we 
do that, the form of the condensate can be simplified further by symmetry arguments. Note that as the Lagrangian 
possesses SU(2) R+L+ c global symmetry and so does the condensate, we should have Ax = A 2 . This is as far as we 
can go before solving for the gaps explicitly using the gap equations. To solve for the gaps we write down the action 
as a function of Ax and A 3 


<S[Ai, A 3 ] 


8 AxA{ + 4A 3 Ag 
4 G 


^Tr(Log (S-^J). 


( 6 ) 


Minimizing the action of Eq. 6 with respect to Ax and A 3 we obtain the gap equations for the two gaps. The gap 
equation for Ax is given by 
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and for A 3 it is 


dS 

dAl 


Ai=A3=A 


4 A 3 
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2 4i I *5quark • 


^quark \ 

dA\ J 


= 0. 
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( 8 ) 


. In order to proceed further we need to evaluate ^quark hr momentum space as the calculations are much simpler in 
it than in position space. However, as is evident from the position space form of S'I u 1 ark given in Eq. 4, the inverse 
propagator appears to not be translationally invariant for the charged quarks, which means that it is not diagonal 
in the momentum space and we lose the advantage we had in going to the momentum space. However, this can be 
remedied by using the Landau level basis for the quarks instead of the standard free particle basis. As the eigen states 
of a fermion in a magnetic field is given by the Landau levels, it is natural to expect the propagator to be diagonal 
in this basis. This method was first introduced by Ritus in the year 1972 [41], We write down the explicit form of 
the quark propagator in the Landau level basis, where l is the Landau level index, in the appendix in order to avoid 
interrupting the thread of logic here. Note that in the appendix we derived the quark propagator only in the limit 
Ai = A 2 = A 3 . However, this will suffice for the purpose of solving Eq. 7 and Eq. 8 as both of them are to be 
solved for Ai = A 2 = A 3 . Also we will be working in the limit of large magnetic fields, of the order of the square 
of the chemical potential so that only the lowest Landau level is occupied. As mentioned before, we simplify 7 and 
8 using the expression for the propagator in the appendix and obtain two equations for A. By demanding the two 
be consistent with each other we can determine the magnetic field as a function of the chemical potential for which 
Ai = A 3 . Simplifying 7 we get 
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and simplifying 8 gives 


1 u 2 , /2A\ 

<? 2 _ vr 2 S (a) 


( 10 ) 


. A is the cut-off used for the divergent integrals in the gap equation. It is worthwhile mentioning at this point that 

although we introduce a sharp ultraviolet cutoff here for the sake of clarity, the cancellation of the UV divergences in 

the calculations that follow do not depend on any kind of regularization procedure. 9 and 10 can only be consistent 

when eB = 2/i 2 . This brings us to the reason why we are considering eB ~ 2/i 2 in the first place and it is because we 

~ 2 

need both high magnetic fileds eB > ^ and all the gaps to be equal to simplify our calculations. The motivation is 
to obtain analytic results in a regime where the effect of the magnetic field on the Goldstone mode propagators would 
be considerable. Let us introduce a matrix at this point given by A = A~| Ai=Aa=A . Having obtained the strength 
of the magnetic field for which the gaps are equal let us now try to find the propagator for the neutral Goldstone 
modes. Before we do that we write them out explicitly here for clarity. Note that the gap transforms under chiral 
flavor rotation as 


_ e *7 5 G a a a ^ e i7 6 (G“) T a a 

= (1 + i-f 5 G a a a ) A (1 + i 7 5 (G a ) T a Q ) + .. (11) 

= A + ia a "/ 5 M a + ..., for a = 1 ,.., 8 

where, G a are generators of SU(3) gauge group, a a are fields corresponding to the generators G a and M a = G a A + 
A [G a ) T . Similarly, under an axial 17(1) rotation which happens to be an approximate symmetry of the Lagrangian 
at high density, the condensate transforms as 

A' = gG 5 ^ 9 A e * 75 “ 9 

(12) 

= A + 2ia 9 7 5 A + ... = A + ia 9 7 5 M 9 + .. 


and under an ordinary 17(1) rotation 


A' = e iaW Ae iaW 

= A + 2ia 10 A + .. = A + ia 10 M 10 + ... 


Among the 10 modes M a we need to find the neutral ones. The neutral goldstone modes are the ones for which 
Q.M a + M a .Q = 0. We find that M 1 , M 2 , M 3 , M 8 , M 9 and M 10 are neutral as expected and the rest are not. We 
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describe the computation of the propagator for a 3 in some detail here and state the result for the rest of the five as 
the calculations are very similar. To proceed we first need to rewrite the action of Eq. 3 in terms of A^ as 
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Now we expand the above action about A = A upto second order in fluctuations(5A) to obtain the following 
expression 
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(15) 


The coefficients of the terms in the above expression which are first order in fluctuations SA go to zero due to the gap 
equation. Also note that there is an easy trick to simplify second derivatives of the trace of logarithm of the inverse 
quark propagator with respect to the gap matrix and it is the following, 


d Tr (Log(5' quark )^ _ ^ ( Q ^quark c ASqu^rk ^ 

d^- jb dA~ kd _ lf - ^{^dA-.^dA^)- 


(16) 


Although we have explicitly written out the above expression for the second derivative with respect to A - , similar 
relations are true for second derivatives with respect to (A - )t and mixed second derivatives with respect to A - and 
(A - )'!'. Now we set SA = 3 in Eq. 15 as we are interested in the propagator of a 3 and denote the part of the 

action that is quadratic in a 3 as £ 3 . After using Eq. 16 and the expressions for the quark propagators for the 0 th 
Landau level from the appendix and some tedious algebra we finally get 
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01 A 


eB 

(17) 


In Eq. 17 we first do the po integral followed by the p and p 3 integrals. The three momentum integrals and 
the p 3 integrals are performed about the Fermi surface, |p| ~ p and |p 3 | ~ p. Note that Eq. 17 has integrals 
which are divergent and we need to carefully subtract these divergences to obtain our final result. The divergences 
basically cancell with the term with the coupling constant (G) in the denominator in Eq. 17. To see how exactly this 
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cancellation works we reorganize Eq. 17 in a form that is more convenient 
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(2tt) 4 2 4 2tt ((p 0 - fco) 2 - (|p 3 - fc 3 | - p) 2 - A 2 ) (p 2 - (|p 3 | - p) 2 - A 2 ) 


Po(Po - fco) - (IP 3 I - M) 2 ” A 2 


dp 3 dp 3 1 47r 2 e-B 16 

"(2^ 2^^2^ «^fc 0 ) 2 - (|p 3 | - p) 2 - A 2 ) (p 2 - (|p 3 | - p) 2 - A 2 ) 


+ ... 


(18) 


. The third and the fourth term in Eq. 18 contain the leading divergences in the expansion Now we need to 
use the gap equation to cancell these divergent integrals with the second term in Eq. 18. Using the gap equation at 
eB ~ 2 p 2 we can see that 


d 4 p 


\eB\ f dp 0 dp 3 1 


4 G 


(2 n) 2 2 (po) 2 — (|p 3 | — p) 2 — A 2 


v (2 tt) 4 (p 0 )2 - (|p| - p) 2 - A 2 4 tt 
which achieves the cancellation we were looking for and we are left with an action for a 3 given by 

d 4 p 


(19) 


<S 3 [A] = <S[A] + A 2 (a 3 ) 


3\2 


3 \ 2 


+ A 2 (a 3 ) 

A 2 (a 3 ) 2 

4 


k 2 

' (2tt) 4 (p 2 - (|p| - p) 2 - A 2 )((p 0 °- fco) 2 - (|p| - p) 2 - a 2 ) 

dp 3 dp 3 l\eB\ kl 

(27T) 2 2 477 (Po ~ (|p 3 | — p) 2 — A 2 )((p 0 — fco) 2 — (|P 3 1 - A*) 2 - A 2 ) 

r d 4 p 


1 + 


P-(P ~ k ) \ 4(po(p 0 - k 0 ) - 4(| p - k| - p)(|p| - p) - 4A 2 ) 

IpI Ip — k| y ((p 0 - fc 0 ) 2 - (|p — k| - p) 2 - a 2 ) (p 2 - (|p| - p) 2 - a 2 ) 


- 1 + 


P-(p-k) 

IpIIp — k| 

A 2 (a 3 ) 2 f d 4 p 
4 


4(po(po ~ k 0 ) - 4(|p| - p)( jpj - p) - 4A 2 ) 

((po - k 0 ) 2 - (|p| - p) 2 - A 2 ) (p 2 - (|p| - p) 2 - A 2 )J 


( 20 ) 


3\2 


- A 2 (a 3 ) 


P-(p-k) 

(277)4 VIpIIp — k| 

dp 0 dp 3 1 |eB| 


- 1 


4(p 0 (Po - k 0 ) - 4(|p| - p) 2 - 4A 2 ) 


((po - ko) 2 - (|p| - p) 2 - A 2 ) (p 2 - (|p| - p) 2 - A 2 ) 
2po(po - ko) - 2(|p 3 - fc 3 | - p)(|p 3 | - p) - 2A 2 


((Po - fco) 2 - (|p 3 ^ fcsl ^ P) 2 - A 2 )(pq - (|p 3 | - p) 2 - A 2 ) 


/ 2(2t7) 2 2 277 
2po(po ^ fco) ^ 2(|p 3 | - p) 2 - 2A 2 


((Po - fco) 2 - (IpsI ^ M) 2 - A 2 )(pq - (|p 3 | - p) 2 - A 2 )J - 

. The above action has only finite integrals except for one term with a factor of ( ^>[|p-k] — l) i ns ide the four 
momentum integral. Although this term is divergent, it is suppressed by powers of ^ and can be ignored for our 
purpose. If we now expand our result in and we obtain the leading order quadratic action in a 3 as 


S 3 [A]=S[A]- 


*((« 3 ) 2 )p 2 

4tt 2 

i((a^)\eB\ 
1677 2 


h 2 _ 
ft. 0 


( 21 ) 


(fco — fc’ 3 ) + 


































































We can perform similar calculations for the other neutral Goldstone modes 1,2,8, 9,10 and find the corresponding 
inverse propagators. From these expressions we can extract the speed of the six neutral Goldstone modes. If we 
denote the speed of the i th excitaion perpendicular and parallel to the magnetic field as v l ± and vl respectively, we 
have = §, (w l )jj = § for * = 1, 2, 3, ( v l )\ = J^, (u*)| = H for i = 8 and (w*)^ = (u*)y = || for i = 9,10 for 
the values of the magnetic field for which eB ~ 2 /i 2 . 


III. LOWER MAGNETIC FIELDS 


The calculation in the section above although applicable to only eB ~ 2/x 2 paves the way for a similar calculation 
in this section for more realistic magnetic fields. In this section we discuss the propagators of the neutral Goldstone 
modes for magnetic fields eB < 0.3/r 2 . Just to clarify the separation of scales in this discussion, if the momentum 
scale at which we are looking at the propagators of the Goldstone modes is denoted by k , we will be concentrating 
in a regime |k | 2 < A 2 < eB < 0.3/r 2 . This regime is such that multiple Landau levels are filled and the total number 


of Landau levels occupied by the quarks is given by 


2 eB 


where 


denotes the greatest integer smaller than 


Before we delve into the details of the effective theory in this regime, it should be noted that the effect of multiple 
Landau levels on the gap equation in the presence of a magnetic field was discussed in [27]. It was found that the 
two gaps Ai and A 3 exhibit de Haas-van Alphen oscillations. This is a consequence of the fact that the number of 
occupied Landau levels changes discontinuously as the magnetic field is varied with respect to the chemical potential. 
From [27] we can see that the two gaps Ai and A 3 become appreciably different from each other for eB > 0.5/x 2 . 
Note that as the magnetic field is increased further, the two gaps become equal to each other at eB ~ 2/i 2 which 
is consistent with our results in the previous section. The effect of the de Haas-van Alphen oscillation is not very 
significant for magnetic fields eB < 0.3/r 2 and it is reasonable to use Ai = A 3 in this regime. As we lower the 
magnetic field further this approximation becomes more accurate[27]. 


Although it is probably difficult to treat this problem in a regime where multiple Landau levels have to be 
taken into account analytically in its full generality, it is possible to do so in the limit |k| <C A < eB, if we 
were to consider only leading order terms in an expansion in -A-. In this limit we can obtain an expression 
for the effective action for the Goldstone modes analogous to 17. The only difference is that in addition to a 


term like — A 


dppdp3 1 47r 2 eB 16 


/ (2tt) 4 2 


Po(Po-fco)-(b 3 - fc 3 |-/d(b3|-A0" A2 


o 


m 


A 2 (a 3 ) 2 
4 


r dpgdp 3 An eB 16 
J (2n) 4 4 2n 2^1=1 


4 2n (( po-fc 0 ) 2 -(|P3-fc3|-M) 2 - A 2 )(pg- ( b3h/r) 2 - A 2 ) 

Po(Po-*o)-(v / b3-fc3l 2 +2|eB|i-Ai) (y/\p 3 \ 2 +2eBl-p) -A 


as in 17, we also have 


°(S) 


where 


(2tt) 4 11 ^(p 0 -fe 0 ) 2 -^ 1v /b3-fe 3 | 2 +2 eBl-pj -A 2 j (p 2 -( 3 /\p 3 \ 2 +2eBl-p,'j -A 2 ^ 

l is the Landau level index. The p 3 integrals in this additional term which corresponds to l 7 ^ 0, are to be performed 

around \Jp1 + 2 eBl ~ /t instead of |p 3 | ~ /t. After subtracting the divergences as before and computing the integrals 
and the sums, we end up with expressions for the speeds of the Goldstone modes along the magnetic field and 
perpendicular to the field. Here we state the results only as the manipulations involved are very similar to the ones 
described in the previous section. The speed of the three modes a 1 , a 2 , a 3 parallel to the magnetic field is given by 


I (eB) = 


V? 1 \eB\ | \eB\ 2y / 2 JiF\ r> p 
12 T 16 ' 16 u 


4 


//, 


_I 1 _ A 

2’ 21 ell | 


+ 


2\eB\ 


\eB\ , \eB\ 2p 
+ 16 v^leA] 


Re 


H) 


Hr 


1 - 


21 ell | 


2|eU| 


- H, 


_I 1 _ A 

2 ’ 21 ell j 


Hr 


1 - 


2|eB| 


( 22 ) 


and 


vj_{eB) = 


4 


til 

12 


|eH| 1 \eB\ 2p p 
16 16 v^jgSf 


(-*) 


H, 


3.1- 


21 ell | 


A 

L 2|eS| J J 


Hr 


I i _ p 2 

2’ 21 ell | 


(23) 


where H^ is the Hurwitz zeta function. Similarly for a 8 we have 

I ell | , \e&\ 2.yJ 2\eB\ 


vl{eB) = 


It¬ 

'd 


48 


48 


-Re 


H 


_I 1 _ 

2’ 21 ell | 


+ 


21 ell | J J 


P 2 

3 


\eB | 
48 


|eB| 2p 

48 


(-0 


H 


1 1 _ p 2 
2’ 2|ell| 


P 2 
[2\eB\ 


-Hr 


_I 1 _ til 

2> 2|eB| 


-Hr 


I 1 _ f 2 

2’ 21 ell | 


(24) 
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0.15 0.20 0.25 0.30 ^2 

(c) 


FIG. 1: Fig. la shows speed in directions parallel and perpendicular to the magnetic fields for the mesonic modes 
1,2,3 as a function of eB where as Fig. lb shows speed in directions parallel and perpendicular to the magnetic 
fields for the 8 th mesonic mode as a function of eB. Fig. lc is a similar plot showing the speeds of mode 9 and 10. 


and 


' (eB) = 


and for a 9 and a 10 


and 


I (eB) = 


hi 

3 


5/x 

18 


5 fi 2 
6 


v 2 ± {eB) = 


5 m 2 
6 


H_ 
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\eB\ . \eB\ 2/t 


48 


+ 


48 




Re 


(-*) 


Hr 


|eB| 

|eB| 2yJl\e.B\ 

3 + 

3 

H 

|eB| 

M| 

2 M 

3 ^ 

3 y/2 \Jb\ 

Mi + 

Ml 

2/2 


Hr 


-u- 


Re 


(-*) 


Hr 


Re 


5/H 

18 


(-9 


ii, 


„2 

D l 

H 2 

2|eB| 

L 2|eB| J 




M 

2|eB| P L 

2|eB|JJ 

.. 2 

, + 

M 2 

2|e£| 

[_ 21 eS | 

a 2 

^ ~ 4- 

n 2 


-Hr 


-Hr 


3.1- 


2\eB\ 


_I 1 __ M* 
2’ 2|eB| 


-Hr 


1 - 


2|eB| 


2|eB| 


II, 


2’ 2|eB|_ 


1 i_M: 

2 > 


2|eB| 


(25) 


(26) 


(27) 


Although the above expresseions are valid for most of 


is restricted 


by (a 


2 eB 


2eH ) > £A 2 and ((1 + 


k| < eB < 0.3/u, 2 , the regime of validity of the above expressions 
gW ^ 2 eB — fj, 2 ') > £A 2 where £ is a number of order one but 


greater than one. We plot the expressions for d_l and dm for all six neutral Goldstone modes as a function of the 
magnetic field in Fig. 1. The empty regions in the curves correspond to a regime where ^/r 2 — 2 eBj < £A 2 


and ((1 + 


2 eB 


2 eB — fj , 2 ) < £A 2 and the above expressions are not valid. In order to determine the regions to be 


excluded in the plot due the reason mentioned above we have taken A ~ 0.1/2. Although it seems from the regions 
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plotted in Fig. 1, that the speed continues to drop to zero as 


2\~eB\ 


2 \eB\ —> p 2 in the regions we have not plotted 


in the figure, this is not the behavior in reality. It is an artifact of how we computed the integrals for the topmost 

r 2 i 

Landau level(Z = lmax = S-fa ) in our calculation. The finite P 3 integrals corresponding to the topmost Landau 


level were performed from + 2\eB\l rnax — n = —oo to y + 2\eB\l rnax — fi = oo which is a reasonable thing to 
do when l 2 — 2eB^j > £A 2 and ((l + ) 2eJ3 — p 2 ^j > £A 2 as the the integrands involved go to zero as 

\Jp\ + 2\eB\l ma x — P > A. In regions (^p 2 — 2eB^j < £A 2 and (^1 + ^ 2eB — p 2 ^j < £A 2 the lower limit 


of integrals for the topmost landau level should be pz = p 2 — 2\eB\l max . The implication of this for the Fig. 1 
seems to be that as the magnetic field is increased, the curve for the speed for l max = l' first decrease as ^g —> l 1 as 


shown in Fig. 1 and then increase over a region of ~ £A 2 about a point where ^g = l', which is not captured in Fig. 
1 , to meet the curve for the speed for l max = l' — 1 . 


IV. CONCLUSION 


In the last two sections we analysed the propagator for the neutral Goldstone modes of the color-flavor locked phase 
in the presence of a magnetic field for both high magnetic fields eB ~ 2p 2 and more realistic moderate magnetic fields 
A 2 < eB < 0.3/j 2 . We found that the speeds of these modes differ considerably from their values in the absence of a 

magnetic field found in the weak coupling limit J |, for the range of the magnetic field we are considering here. There 


are a couple of improvements that can be made to our calculations here. The first one is for the regime eB ~ 2p 2 . The 
magnetic field in this regime is so high that only the lowest landau level is filled. This may imply that the magnetic 
field should already be probing the vacuum which can lead to the condensation of the chiral condensate ( qq ). This 
phenomenon at zero density is known as ‘magnetic catalysis’ [42-47],[48-52],[53, 54]. This is likely to have an impact 
on the magnitude of the MCFL gap. If this is the case, then the mesonic excitations will involve quark-antiquark 
quantum numbers like mesons at zero density along with the diquark quantum numbers considered here. This will in 
turn add corrections to the anisotropy in the speed of the neutral mesonic excitations that we found here. Also our 
calculation at eB ~ 2 p 2 involves gaps Ai = A3 = A which is why although we can accurately (barring the corrections 
coming from the presence of a chiral condensate) predict the anisotropy in the speed of the Goldstone modes at 
eB = 2p 2 , we lose information regarding how this anisotropy changes as a function of around the point eB = 2p 2 . 
In order to recover this information the calculation for the propagators of the Goldstone modes needs to be carried 
out in its full generality using gaps that are unequal or Ai ^ A 3 . This is not very important at lower magnetic fields 
eB < 0.3 p 2 as the two gaps are close to each other in magnitude in that regime. However, for large magnetic fields 
2eB > p 2 this is important as the gaps become vastly different [27] as the magnetic field is increased from eB ~ ^ 
onwards, except for the point eB ~ 2 p 2 . A more detailed calculation including the effect of magnetic catalysis and 
with unequal gaps Ai and A 3 will improve the accuracy of this calculation and will also give us a prediction of how 
the anisotropy in the speed of the Goldstone modes changes as the magnetic field is varied in this regime of very high 
magnetic field. 

As seen in Fig. 1 the speed of the Goldstone modes not only differs considerably from the standard weak coupling 
value of that is found in the absence of a magnetic field but also undergoes de Haas-van Alphen oscillations with 

2 

increasing magnetic field for fields eB < These could have observable consequences for the transport properties 
of the core of a magnetar. The biggest contribution to the neutrino emissivity of the CFL phase comes from the 
decay of the pseudo Goldstone modes or the scattering of the Goldstone modes into neutrinos. The anisotropy found 
in the speed of the Goldstone inodes here could significantly alter the neutrino emissivity parallel to the magnetic 
field compared to that perpendicular to it. Also the oscillations observed in the speed could inroduce oscillatory 
behavior in the neutrino emissivity which could be relevant for a magnetar going through stellar evolution. It was 
found in [27] that during its stellar evolution, a magnetar with a CFL core could go through several successive phase 
transitions which would induce discontinuous changes in the induced magnetic field in the core. This could translate 
into observable effects in the cooling due to oscillatory nature of the neutrino emission. 
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VI. APPENDIX 


The Landau level basis has been explained in great detail in [39, 55]. Here we outline the basic steps to find the 
quark propagator for the problem at hand. We only need to find the propagator in the limit Ai = A 3 . In its current 
form neither A - nor its square is diagonal. However our calculations become a lot simpler if the square of the gap 
matrix is diagonal in the basis we are working in. After we find such a basis for A~ | A A the condensate in this 
new basis is given by 


I Ai=A 3 


where V is the matrix given by 


and A' is, 


A' = 



A' 

= V 

(A- 

Ai = 

As) 


/ / 

2 0 

3 u 

0 0 

1 

76 

0 0 


0 

1 

0 0 

0 

0 0 


0 

0 

1 0 

0 

0 0 


0 

0 

0 1 

0 

0 0 

v = 

0 

0 

0 0 

1 

72 

0 0 


0 

0 

0 0 

0 

1 0 


0 
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0 0 

0 

0 1 
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0 

0 0 

0 

0 0 


1 

V 73 

0 

0 0 

1 

75 

0 0 

/-A 

0 

0 

0 

0 

0 

0 

0 

0 

-A 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-A 

0 

0 

0 

0 

0 

0 

0 

0 

-A 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-A 

0 

0 

0 

0 

0 

0 

0 

0 

-A 

V 0 

0 

0 

0 

0 

0 


V 


(28) 


J_ \ 

7e 

0 

0 

0 

1 

\/2 

0 

0 

0 

75 ) 


(29) 


0 

0 

-A 

0 

0 

0 

0 

0 

0 


0 \ 
0 
0 
0 
0 


-A 0 
0 0 
0 0 
0 2 A/ 


(30) 


with Ai = A 3 = A. We need to evaluate the quark propagator in order to proceed. The quark propagator and its 
inverse are two by two matrices in the Gorkove space, 


S, 


quark 


/ oil C 12 \ 

_I '"'quark '■'quark | 

“ I o21 o22 I 

y '' 7 quark ^quark'/ 


(31) 


'S’quark’ *^quark> ^qiarki ^quark are all 9 x 9 matrix in the color flavor space. The inverse quark propagator ( S quar kY 
is diagonal in the Landau level basis and can be expressed as, 


((Squark) -1 ) = ((p° + A*b° - Wrf - p 3 7 3 )D+ + ((p° + /i) 7 ° + - p 3 7 3 )D_ + ((p° + /z) 7 ° - p 3 7 3 l 

> 

((•S’quark) -1 ) = ((p° - m)7° - V2eBlj 2 - p 3 7 3 )D_ + ((p° - /z) 7 ° + V2 eBlj 2 - p 3 7 3 )D + + ((p° - p)j° - _p 3 7 3 ()S5) 
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((S^k)- 1 ) 12 =i'f*' = K (34) 

and 

((^uark)” 1 ) 21 = 7°(-*7 5 A')7° = K- (35) 

We have defined two matrices A' c and A' c as we need them for convenience later in the appendix. We need to invert 
‘-’quark to § e t the quark propagator Sq Uar k in the landau level basis. The inverse of the elements of 5 quar k are given by 


(SqL-k) 1 


((p° + m)7 0 - VFFlj 2 - p 3 7 3 )^+ + ((p° + p)l° + VFFlj 2 - p 3 "/ 3 )tt^ + ((p° + m)7 ° - p 3 7 3 )^0 
- (? 7 5 A') \{{p° - fi)j° - V2£Bl 1 2 -p 3 7 3 )fl_ + ((p° - ^ 7 ° + \[FFl 1 2 - p 3 j 3 )tt + 


+((P° - Mb° - P 3 7 3 )^o] 1 (7°(-i7 6 ) (A') f 7°) 


(36) 


(Squark )- 1 


((p° 


p) 7 ° — V 2 eBlj 2 
(7°(-i7 5 )A ,t 7°) 


- p 3 7 3 )f2_ + ((p° - p)7° + V / 2d^7 2 - P 3 7 3 )^+ + ((p° - p)7° - P 3 7 3 )^o 
[((p° + ph° - VFFl'y 2 - p 3 7 3 )fl+ + ((p° + Ai)7° + VFFl'y 2 - p 3 7 3 )f2_ 


+((P° + M)7° - P 3 7 3 )^o] 1 (*7® A') 


(37) 


(^quark )” 1 


('C 21 'l-l 
V^quark/ 


((p° + m)7 ° - V / 2gIi 7 2 - p 3 7 3 )^+ + ((p° + M)7° + FFFlj 2 - p 3 7 3 )fl_ + ((p° + M ) 7 ° - p 3 7 3 )^ 0 
(* 7 5 A') 5 22 ]" 1 

(38) 


((p° - m)7 ° - F~2Wrf - p 3 7 3 )fl_ + ((p° - m)7 ° + VFFrf - p 3 'y 3 )fl + + ((p° - M ) 7 ° - p 3 7 3 )flo 
(7°(-i7 6 )A / t7°)5 11 ]- 1 . 

(39) 


To write out the different components of S quar k explicitly we need to define the following projectors in the massless 
quarks limit, 


1 / ± 7 °(7 3 p 3 + 7 2 \/2iBl) \ 

2 V sjpl + ^Bl ) 

1 f {± 7°(7 : V-7 2 V / 2eB/)\ 

2 V sjpl + TeBl ) 


(40) 


(41) 


From this point onwards we refer to S'quark as S. We write the diagonal terms of the propagator in the Gorkov spacep 
in terms of these projectors as 


/Qll\ _ oil _ oil _ oil 

\ D Jll-ll — J 12-12 ^ d 21-21 — J 22-22 


(p 0 - |pI - m)7°A + (pq + jp| - A<)7°A 

p 2 - (|p| + p) 2 - A 2 + pg - ( |p| _ m) 2 _ A 2’ 


(42) 
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nil _ oil _ 
°13-13 — °23-23 — 


oil _ nil _ 
°31-31 — °32-32 — 


(Po - \/pi + ZeBl - p) 7 °At ^ (Po + \Jpi + 2 eBl - p) 7 °A_ 

Po - ( \JpI + 2 eBl + n) 2 - A 2 Po ~{\]pI + 2 ^l - p) 2 - A 2 

(Po ^ 'Jpl + ZeBl - p) 7 °A+ ^ (p 0 + Y^pi + 2 eBl - p) 7 °A“ 
Po - (\/p§ + + A ) 2 - A 2 Po - (V // Pi+ _ 2A^“ p ) 2 - A 2 


(43) 


(44) 


and 


b 1 — 

3 33-33 ~ 2 


(po - |p| - p) 7 °A + , (po + IpI - p) 7 °A _ 


Po - (IpI + p ) 2 -4A 2 pg- (|p| -b 2 -4A 2 ’ 


?22 


qzz _ qzz _ nZZ _ n. _ 

° 11-11 — ° 12-12 — ° 21-21 — ° 22-22 — ~ 2 ~ 

Po 


(po - IpI + p) 7 °a + , (p 0 + |p| + pb°A _ 


-p) 2 -A z p 2 -(|p|+p) 2 -A^ 


s 


-22 _ o22 _ 

13-13 — °23-23 — 


o22 _ n'. 

Oqi_qi — 0 ‘ 


_ ,22 _ 

'31-31 — °32-32 


(p 0 - \JpI + 2 .eBl + p) 7 °A+ ^ (po + \Jp\ + 2 eBl + p)^°A + 

3 o - (\jpi + 2 eBl - P ) 2 ” A 2 Po ~{\JpI + 2 ABZ + p ) 2 - A 2 

— pi + 2 eBl + p) 7 °A+ (p 0 + \jpi + 2eBl + p)"/°AZ 

/ — I / = ’ 


(Po 


pi - {^f+2iFl - p ) 2 - A 2 Po - (\/p§ + 2ePZ + p ) 2 - A 2 


(45) 


(46) 


(47) 


(48) 


q22 

°33-33 — 2 

Po 


(po - IpI + pb°A + , (p 0 + |p| + pb°A~ 


— fi ) 2 — 4A 2 Po — (IpI + p ) 2 — 4A 2 


(49) 


where the upper indices on 5' are Gorkov indices and lower ones are color-flavor indices. To write down the expressions 
for S ' 12 and S 21 with some clarity we need to define the following 


7 22\-l _ 7°A+ 


{Z zz )~ = 


7°A~ 


Po a IpI p po — |p| — pt ’ 


(50) 


(M ZZ )~ L = 


7° At 


7°a: 


Po + y 2eBl + pg — p po — y 2ePZ + P 3 — A 1 


(51) 


(P 22 )"i = 


7° At 


7° A] 


Po + y 2ePZ + P 3 — p Po t/ 2ePZ + p§ — p 


(52) 


(b 1 )" 1 EE 


7°A+ 


7°A- 


po + IpI + m Po — IpI + m ’ 


(53) 


(M LL y = 


7° At 


7°a: 


Po + y 2eBl + p| + p Po — y 2ePZ + P 3 + A 1 


(54) 


(P 11 )" 1 ee 


7°A+ 


7 °a: 


Po + y 2eBl + pi+ p po — y 2ePZ + P 3 + A* 


(55) 
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S 21 and 5 12 can be written in terms of these newly defined quantities as 


: -((S- 1 ) 22 ) _1 (7°(-i7 5 )A / t T 0 )S 11 
-(z 22 )~ 1 a' c sH 0 


~(Z zz )~ L A' C S{ 
0 


0 

-(Z 22 )- 1 A' c sl 
0 


-(M zz )~ l A' c S{ 


-(z zz )~ 1 a , c s£ 


_(P ZZ ) — L A' c Sk 


(56) 


-(M 22 )- 1 A' c S\l 
0 


-{p 22 )~ 1 a , c sI\ 

0 


2(Z 22 )- 1 A' C S\\) 
(57) 


and 


f-(Z 11 )~ 1 A , c S 22 
0 


-(Z LL )~ L A' c Sf 


-(z ll )- l a' c s$ 


~(P 11 )- 1 A'S 2 


-(Z 11 ) -1 A' C S% 
0 


-(M 


-(P 1 


) _± A' C S| 

0 


' c s 22 

0 


(58) 


-(M 11 )~ 1 A' 


2(Z 11 )~ 1 A' c S%p 

(59) 


where S stands for . 

L J L J L J 
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